Distribution of continuous-variable entanglement by separable Gaussian states 
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Entangling two systems at distant locations using a separable mediating ancilla is a counterintu- 
itive phenomenon proposed for qubits by T. Cubitt et al. [Phys. Rev. Lett. 91, 037902 (2003)]. 
We show that such entanglement distribution is possible with Gaussian states, using a certain 
three-mode fully separable mixed Gaussian state and linear optics elements readily available in ex- 
periments. Two modes of the state become entangled by sequentially mixing them on two beam 
splitters, while the third one remains separable in all stages of the protocol. 



Quantum entanglement is a striking property of com- 
posite quantum systems that lies at the heart of the fun- 
damental quantum information protocols such as quan- 
tum teleportation [l[ or quantum cryptography [2] . The 
typical scenario involves two parties in distant laborato- 
ries, Alice holding the quantum system a and Bob in pos- 
session of the system b. They need to establish a quan- 
tum channel between their remote location in the form 
of the shared entangled state, which cannot be prepared 
merely by local operations on systems a and b and clas- 
sical communication between Alice and Bob Having 
only separable systems at hand and not having a pos- 
sibility to meet each other in one place, Alice and Bob 
can entangle their distant quantum systems only by em- 
ploying another ancillary quantum system c, which first 
couples with the system a, then is send over to the re- 
mote location where it interacts with b. This is a required 
global operation that facilitates entanglement between a 
and b. Remarkably, a can be entangled with b by send- 
ing the ancilla c that becomes never entangled with the 
subsystem (ab) Q . In [i[ the counterintuitive effect of en- 
tanglement distribution by separable ancilla was studied 
in the context of finite-dimensional systems. 

In this paper we show how to turn this idea into a prac- 
tical concept. We consider infinite-dimensional quantum 
systems, e.g., light modes. We propose a feasible three- 
step protocol where two distant separable modes A and 
B become entangled after interacting stepwise with the 
third mode C. At any stage of the protocol, the mode C 
is separable from the subsystem (AB). Our scheme relies 
entirely on Gaussian states and the challenging nonlin- 
ear controlled-NOT gates of the previous idea [J] are re- 
placed by simple beam splitters. Therefore, the protocol 
can be implemented with Gaussian states and operations 
that are currently available in the laboratory. Moreover, 
the proposed protocol allows a more simple determinis- 
tic distribution of entanglement than the previous qubit 
protocol that requires an additional operation on systems 
b and c on Bob's side 

Our protocol is schematically depicted in Fig. [TJ The 
aim of the protocol is to entangle mode A in Alice's lab- 
oratory with separable mode B in Bob's distant labora- 
tory by sending a separable mediating ancillary mode C 
from Alice to Bob. For pure quantum states this is not 
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FIG. 1: Scheme of the protocol for distribution of continuous- 
variable entanglement by separable Gaussian states. Step 1: 
LOCC preparation of a fully separable Gaussian state of three 
modes A, B and C. Step 2 entangles mode A with a pair 
of modes (BC). Step 3 entangles mode A with mode B. 
Mode C remains separable from the pair of modes (AB) in all 
steps. BSac.bc denote balanced beam splitters, Da,Db, Dc 
are local displacements distributed according to the Gaussian 
distribution with correlation matrix Q. 



possible [4J. Therefore, Alice and Bob have to construct 
by local operations and classical communication (LOCC) 
a suitable mixed fully separable Gaussian state of three 
modes A, B and C. The engineering of such a state is the 
aim of step 1 and represents the most challenging part of 
the problem. Alice and Bob start with three pure single- 
mode Gaussian states. They prepare modes A and B in 
the same momentum-squeezed vacuum states and rotate 
them clockwise and anticlockwise, respectively, by the 
same suitable angle. The ancillary mode C is initially on 
Alice's side and is in a vacuum state. Alice and Bob then 
displace locally the three modes by random correlated 
displacements Da,Db and Dc with Gaussian distribu- 
tion characterized by a correlation matrix Q(x) specified 
below. This procedure generates the desired three-mode 
mixed Gaussian state with separability properties tun- 
able by changing the parameter x. 

The actual entanglement distribution commence in 
step 2. By mixing modes A and C on a balanced beam 
splitter BSac Alice entangles mode A with the pair of 
modes (BC) while mode B is separable from (AC) and 
mode C is separable from (AB). Next, she sends mode 
C to Bob. In step 3 Bob mixes modes B and C on a 
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balanced beam splitter BSbc finally entangling A and 
B, while C still remains separable from (AB). 

The modes A, B, and C are described by three pairs 
of canonically conjugate quadrature operators Xj,Pj, j — 
A, B, C. The operators satisfy the canonical commuta- 
tion rules that can be compactly expressed as = 
-iCljk, where £ = (xa,Pa,x b ,Pb,x c ,Pc) T is the vector 
of quadratures and fl = ©f =1 J is the symplectic matrix, 
where J = — ia y (a y denotes the y Pauli matrix). Quan- 
tum states of three-mode system can be represented in 
phase space by the Wigner function [j| of six real vari- 
ables and Gaussian states are defined as those having 
a Gaussian-shaped Wigner function. Any three-mode 
Gaussian state p is therefore fully characterized by the 
vector of first moments £ = Tr(p£), that we assume to be 
zero, and by the 6x6 real symmetric covariance matrix 
(CM) 7 with elements j jk = Tr (p{^ - ^1, & - £ fe l}), 
j,k = l,...,6, where {A, B} = AB + BA. 

Preparation of the three-mode fully separable state. 
We start with the three-mode Gaussian state, which CM 
is composed from the CM of an entangled state jab 
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(1) 



and a noise term in the form of a nonnegative multiple 
of a positive semidefinite matrix P = q\q[ + <Z29j : 



71 (x) = jab ® 1c + x(qiqf + qiq^), 



(2) 



where x > 0. The parameters involved in the CM ([TJ) 
are given by a — cosh(2r), c = sinh(2r) and we assume 
d > r > 0. This is a two-mode squeezed vacuum state 
with the squeezing parameter r with modes A and B 
squeezed, in addition, by local squeezing operations Sa — 
Sb = diag(e d , e~ d ). Our design of the noise term in 
Eq. ((21) is inspired by the method Q used to construct 
various three-mode entangled Gaussian states and q\ , qi 
read 



qi = (O,shi(/.,O,-sin0,\/2,\/2) T , 
q 2 = (cos</>, 0, cos</>, 0, V2, \/2) T , 



(3) 



tan<£ = e Zr sinh(2d) + <J 1 



sinh 2 (2d) 



with sin <f>, cos 4> > 0. This additional noise is chosen such 
that for sufficiently large x the CM ((2]) describes a fully 
separable state. 

The state described by CM |T]) can be naturally pre- 
pared by mixing on a balanced beam splitter Uab |2l 
modes A and B, each in a pure momentum-squeezed vac- 
uum state with the variances of the position quadratures 
({Ax A ) 2 ) = e 2{ - d -^ and <(Ax B ) 2 ) = e 2 ( rf+r ) respectively. 
The entire three-mode state with CM ([2]) then can be cre- 
ated by adding a vacuum mode C with CM tc to the CM 
jab and performing local random correlated displace- 
ments of modes A, B and C distributed with Gaussian 



distribution with correlation matrix xP Q. Making use 
of the criterion of full separability for three-mode Gaus- 
sian states [6j] one then finds that for all x > x scp , where 



2 sinh(2r) 



(4) 



where 5 = e 2d sin cf) + e~ 2d cos 2 <f>, the CM ^ describes a 
fully separable state. However, although the state is fully 
separable the way of its preparation described above is 
not suitable for our purposes. Namely, it is not prepared 
by LOCC but instead requires Alice and Bob to meet 
to implement the beam splitting operation Uab on their 
modes A and B. 

Still the CM ((2|) corresponds to a fully separable state 
and therefore there exists a recipe how to create this state 
by LOCC. The recipe is based on the three-mode sepa- 
rability criterion 0, Q according to which a three-mode 
Gaussian state with CM (x) is fully separable iff there 
exist single- mode CMs 7^, Jb and 7c such that 

Q(x)= 1i {x)-~/a®Jb®Jc>0- (5) 

Interestingly, such single-mode CMs can be indeed found 
for x > x scp in the form 



where 



1A,B 



-2r 



a + (3 =Ft 
=Fr a - 13 



7c = 1, 



(6) 



[e 4r + cosh(4d) - sinh(4<f) cos(2<£)] , 
(3 = { [e 4r - cosh(4d)] cos^) + sinh(4d)} , 

sin(20), (7) 
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and the parameters satisfy the purity condition a 2 = 
(3 2 + t 2 + 1. The CM 7 C represents a vacuum state. 
The CM 7^ (j B ) corresponds to the pure momentum- 
squeezed vacuum state with squeezing parameter s = 
4 In (a + Vet 2 — l) rotated clockwise (anticlockwise) by 

the phase = arctan 

It remains to show that the matrix Q(x) is positive 
semidefinite for x > a; SC p. It is sufficient to show that for 
x = x scp since if Q(x scp ) > 0, then Q{x) = Q(x sep ) + (a 



C p)P > for all x > x scp because (x 



,)P is also 



positive semidefinite. To obtain the eigenvalues of the 
matrix Q(a;sop), we will calculate the eigenvalues of the 
matrix UabQ(x SC p)U^ b , which possesses the same eigen- 
values. They read explicitly as Ai^.3,4 = 0, A5 = 9x sep 



and A 6 = (e sin + e~ 4<1 cos cp) x scp . All of the eigen- 
values are nonnegative and therefore the matrix Q{x) for 
x > a; SC p is indeed positive semidefinite. 

Creation of the fully separable state with CM 71 (x), 
where x > x sop , is now straightforward Q. Initially, Alice 
prepares in her laboratory mode A in a pure single-mode 



3 



squeezed state with CM ^a and the ancillary mode C 
in the vacuum state. Similarly, Bob prepares the mode 
B in a pure single- mode squeezed state with CM jb- 
In the next step, Alice and Bob displace locally their 
modes by random correlated displacements distributed 
according to the Gaussian distribution with correlation 
matrix Q(x). As a result, they prepare by LOCC a three- 
mode fully separable Gaussian state with CM @. For 
the sake of simplicity here and in what follows we do not 
write explicitly the dependence of CMs on the parameter 
x and we implicitly assume that x > x scp . 

Entanglement distribution. In step 2 Alice superim- 
poses modes A and C of a fully separable state described 
by the CM 71 on a balanced beam splitter Uac that 
transforms the CM as 



72 = UacIiUac- 



(8) 



Apparently, the CM is separable with respect to par- 
tition B — (AC). More interestingly, mode C can re- 
main separable from the subsystem (AB) if we choose 
the parameters d,r and x properly. To prove this, we 
apply to CM JHJl the separability criterion based on the 
symplectic invariants [9j. The criterion utilizes the ma- 
trix 7^ c = Ac7 2 A c , where A c = diag(l, 1, 1, 1, 1, -1), 
that describes CM 72 afterpartial transposition with re- 
spect to the mode C 0, S, M, ED, [ij, [13]. The ma- 
trix "f2° nas three symplectic invariants denoted I\ , I2 
and I 3 = det(7 2 ) that can be obtained as coefficients 
of the characteristic polynomial of the matrix ^7^°, i.e. 
det(^7 2 rc -yl) = y 6 + l x y A + I 2 y 2 + 1 3 . According to the 
criterion mode C is separable from modes (AB) iff 



h 



h - 1 > 



(9) 



holds [l4j ■ Calculating the invariants using the equations 
above we arrive at the following simple expression: 



£ = x(ux + v), 



(10) 



where u and v are complex functions of parameters d and 
r given elsewhere (l5j . In this paper we are interested 
in demonstrating the possibility of entangling A and B, 
while keeping C separable. So we show that for particular 
values of parameters d and r we get £ > and hence CM 
72 is separable for x larger than a certain threshold value 
x t h- Eq. pop determines a parabola in the (x, £) plane 
that intersects the x-axis in the origin. Taking e 2 ( d ~ r ) = 
3/2 and e 2< - d+r ^ = 2 (which corresponds to e" 2s w 0.6387 
and 9 « 5.73°) we get u > and v < and the parabola 
is oriented upwards. Making use of Eqs. Q and (TIT)]) , 
the threshold value x t h — —v/u rss 1.04 > x scp « 0.2043, 
and hence for 2 > xth we have E > and the CM (JHJ) is 
separable with respect to the partition C — (AB). 

The protocol is finalized in the step 3. After receiving 
mode C from Alice, Bob superimposes this mode with 
his mode B on another balanced beam splitter BSbc 
■ The CM of the resulting state reads 



73 = U B c72Ub C . 



Remarkably, 73 exhibits entanglement between modes A 
and B whereas mode C remains separable from (AB). 

To verify entanglement between modes A and B we 
express the two-mode CM 73,^5 of the reduced state of 
modes A and B in the block form 



13,AB 



A C 
C T B 



(12) 



with the submatrices A, B, and C of the form: 



r ( c + + (9ihi - -j^)x 



■(ho 



7^ 



c- - (g h Q + -fyx J ' 



A = 
B = 



a + + (di + l)^ (.90 + 9\)x 
(g +gi)x a- + (g 2 + l)x 

b+ + (h\ + \)x 



V2 



V2 

(hi 



(13) 



±2d r 



l)/2;6± = [e ±2d (3aT2\/2 C ) + i]/4; 



where a± 

c± = [e ±2d (a=F\/2c)-l]/2V2; g, = 1 + sim> + jf )/v^; 

= [V2-(-l)J] sin(0+i|)/2+(-l)VV2, j = 0, 1. The 
entanglement of CM 73^5 can be proved if we calculate 
the so called symplectic eigenvalues [l6[ of the matrix 
tIab = A 2,s73,abA 2 ,b, where A 2 ,b = diag(l, 1, 1, -1) 
17| . The matrix has two symplectic eigenvalues v, v' 
that can be computed from the eigenvalues of the matrix 
^2j3AB' wh ere ^2 = ©f=i J and are equal to ±iv'} 
[l7T | . The mode A is entangled with the mode B iff v < 1 
or v' < 1 . The lower symplectic eigenvalue of the matrix 
Ti^as can be- expressed as [ljj 



K- a/k 2 - 4det(7 3ij4 s) 



(14) 



where k = det(A) + det(B) - 2det(C) and ^4., 6, C are 
defined in Eq. (p~3|) . By the same token as in step 2, we 
revert to the case e 2{ ~ d - r ) = 3/2, e 2 ( <i + r ) = 2. Taking 
a; = 1.041 > xth, we obtain an exact expression for v in 
terms of square roots that approximately equals to v w 
0.9571 < 1, which is a clear evidence of the entanglement 
between modes A and B. Note, that by measuring a 
quadrature (xc+Pc)/ V% on mode C the eigenvalue can 
be further reduced (i.e. entanglement can be increased 
[3) to v m = 0.9421. 

Finally, we have to show that the ancillary mode C is 
separable from the two- mode subsystem (AB). We use 
again the simplectic invariants criterion of [9j . Analogous 
to step 2, we calculate the characteristic polynomial of 
the matrix ^7^° and find three symplectic invariants Ji , 
J 2 and J 3 = det(73), which have to obey the condition 
E = Ji - J2 + J 3 - 1 = x(wx + z) > (c.f. ©, (HOD). 
The parameters w and z are again complex functions 
of the parameters d, r [15| . Assuming e 2 ( d ~ r ^ = 3/2, 



?2(d+r) 



2 and x = 1.041, we obtain E « 0.3957 > 0. 



(11) Thus the mode C is separable from (AB). 
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FIG. 2: Performance of our protocol in dependence of the 
variances <(A:ta) 2 } = e 2( - d ~ r) and {{Axb) 2 } = e 2(<i+r) . Gray- 
scale region: all assumptions of the protocol are satisfied and 
v < 1. White region to the right of the slash: either v > 1 or 
C is entangled with (AB) in some stage of the protocol. The 
contour lines display the values of symplectic eigenvalue v. 

In an experiment, verification of entanglement of 
modes A and B in step 3 can be done by measuring the 
entire CM 73,ab similarly as in and applying Simon's 
[l2| or Duan's [l3| separability criterion. The separability 
of mode C from the pair of modes (AB) in steps 2 and 3 
can be proved by measuring the three-mode CMs 72 and 
73. For Gaussian states one then can use the positive 
partial transposition criterion (kJ EH that is sufficient 
for separability of these 1 x 2-mode systems Q. 



As a numerical evidence of the robustness of the proto- 
col, we obtained eigenvalue v — 0.9787 for 7+ + 2 x 10~ 2 1 
corresponding to the initial CM perturbed by a weak 
isotropic noise. Numerical analysis also verifies the per- 
formance of the protocol for a broad range of variances 
((Axa.b) 2 ) = e 2 ( dT '') depicted by a gray region in Fig. [2] 

In conclusion, we have demonstrated the possibility to 
distribute entanglement without sending entanglement 
in infinite-dimensional systems. Remarkably, one can 
entangle two distant modes by a separable mode us- 
ing experimentally feasible Gaussian states and opera- 
tions involving single-mode squeezed states, correlated 
displacements and beam splitters, dispensing with the 
CNOT gates of the qubit case. The distributed entan- 
glement is distillable [2(3] and therefore can be used for 
quantum communication. In contrast with two qubits, 
two light modes can be entangled deterministically even 
without any additional operation on Bob's system and 
ancilla beyond step 3 (cf. 4]). Furthermore, we have 
elaborated the procedure to design three-mode Gaus- 
sian states with desired separability and noise proper- 
ties. Together with the possibility to distribute distil- 
lable continuous-variable entanglement without sending 
it through the channel, it prepares the ground for bet- 
ter understanding and engineering of optical quantum 
networks, continuous-variable cryptography and other 
entanglement-based communication protocols using light 
modes and/or atomic ensembles. The support of the EU 
project COVAQIAL (FP6-511004) under STREP and 
the Czech Ministry of Education (Grant Nos. MSM 
6198959213 and LC06007) is acknowledged. 
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